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Abstract 

We present a new purely equilibrium microscopic approach to the descrip- 
tion of liquid-glass transition in terms of space symmetry breaking of three- 
and four- particle distribution functions in the cases of two and three dimen- 
sions, respectively. The approach has some features of spin glass theories as 
well as of density-functional theories of freezing. 

The main purpose of the report is to present a new purely equilibrium microscopic 
approach to the description of liquid-glass transition in terms of space symmetry 
breaking of three- and four- particle distribution functions in the cases of two and 
three dimensions, respectively. The approach has some features of the spin glass 
theories as well as of the density-functional theories (DFT) of freezing. 

It is usually believed that there are two essential differencies between spin glasses 
and real structural glasses: 1) in the Hamiltonian of spin glasses there is explicit 
randomness from the very beginning, while in the case of real glasses there is no such 
randomness. 2) In experiments with spin glasses there is always the range of the con- 
centration of magnetic impurities where nothing else that a spin glass phase appears 
while in the case of space glass there exists a crystalline ground state. However, in 
real systems one can consider these differencies simply as time scales differencies for 
the freezing of corresponding degrees of freedom with respect to the time scale of 
the real or computer experiments. In fact, there are now some indications that two 
possible candidates for equilibrium glasses do exist: some polydisperce hard-sphere 
systems and some binary mixtures of hard spheres. Even if it is not so, it seems 
us that one needs an "underlying" equilibrium theory of liquid-glass transition to 
understand what really glasses present as space symmetry breaking problem. We 
should mention that beautiful and fruitful time-dependent mode-coupling theory 
[I]] which describes a number of subtle experimental facts does not consider the 
problem of space symmetry breaking. Some other arguments can be found in the 
recent papers by Parisi (see, e.g. and references therein). 

To describe different kinds of space symmetry breaking we use the formalism of 
classical many particle conditional distribution functions 
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Here F s (rx, r s ) is the s-particle distribution function. The functions F s+ i(ri|r ( j ) ...r°) 
satisfy the equation |§ 
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Here z is the activity, p is the mean number density, Sfc + i(ri, r^+i) is the irre- 
ducible cluster sum of Mayer functions connecting (at least doubly) k + 1 particles, 
(3 = l//cgT and T is the temperature. 

If one takes the derivative of (P relative to r 1; one obtains the equilibrium 
Bogolubov hierarchy [[| along with the explicit expression for F s+ 2 as the functional 
on F s+ i which gives the formally exact closure. However it contains infinite series 
and integrals and one has to use some approximations to exploit it. The same can 
be said about the Eq.(|]) itself. 

Let us now consider the symmetry breaking of the one-particle distribution func- 
tion and formulate briefly DFT of freezing (see || and the reviews 0). The equa- 
tion (p for s = is the extremum condition for the free energy functional of the 
inhomogeneous system with the density p(r) = pF\(r) and has the form: 
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or 



Tlk B T = J rfr 1 p(r 1 )[ln(A d p(r 1 ) - 1] - ^.^(r)]/^. (3) 

The excess free energy jF ex [p(r)]/fcsT is just the generating functional for direct 
correlation functions 

c „ (r ,.. r „) = (4) 
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so that Taylor expansion around the liquid can be written in the following form: 
(3AF = J drg(r)\n^ -Y^yj c (n) (r 1 ,...,r k )Ag(r 1 )...A Q (r k )dT 1 ...dv k , (5) 
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where 

Ag(r) = g(r) - g t 

is the local density difference between solid and liquid phases. 

The full system of equations to be solved in DFT contains the nonlinear integral 
equation for the function p(r), obtained as the extremum condition for the free 
energy and the equilibrium conditions for the chemical potential and the pressure 
written in terms of the same functions as in (D). To proceed constructively in the 
frame of DFT we must choose a concrete form of the free energy functional - a 
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kind of closure or truncating - and we must make an ansatz for the average density 
of the crystal. The importance of such an ansatz follows from the fact that we 
are dealing with a theory which is equivalent to Gibbs distribution and one has to 
break symmetry following the Bogoliubov concept of quasiaverages ||. Now it is 
necessary to specify the crystal symmetry (e.g. lattice type) and to locate the freezing 
transition for that particular lattice type 

Ap(r) = pi E k <A<e ikr = pm + Pi<p(t), , r . 

The sum is over reciprocal lattice vectors and the integral is taken over the elemen- 
tary lattice cell A .(p^ are the order parameters of the problem. The DFT approach 
occurs to be very fruitful and was used to calculate a lot of melting curves for 
different systems. 

The 3D DFT scenario of freezing is valid for some 2D systems. However, there 
is a number of 2D systems which melts through two continuous phase transition 
including intermedeate (so called hexatic) anisotropic liquid phase. The scenario for 
such a case of 2D melting is the well known KTNHY || phenomenological scenario. 
We develop a microscopic approach to 2D melting || [LCj in the spirit of 3D DFT. 
Our approach differs from the standard DFT theory of freezing in two main points: 
First, we allow the Fourier coefficients Pg{ t ) of the one-particle distribution function 
expanded in a Fourier series in reciprocal-lattice vectors {G}: p(r) = J2g PG( r )e iGr 
to fluctuate and to have amplitude and phase. Second, we allow the liquid to be 
anisotropic: we consider as possible the existence of a phase with constant density 
but angular dependent two-particle distribution function F 2 (r 1 — r ) 7^ g(r 1Q ). 

These two points of generalization define the two new order parameters: the 
fluctuating Pg{ t ) an d the Fourier coefficients characteristic for the broken symmetry 
of the function F 2 (ri — r ). Our approach again is based on the Eq.(|I]) but now, 
considering hexatic phase, we are dealing with the bifurcation of the solution for 
the two-particle distribution function. The relative spatial distribution of pairs of 
particles is characterized by the function F 2 (ri\r ) = F 2 (ri — r ). The vector r x — r 
defines the direction of the bond between the molecules at the points ri and r . In 
the ordinary isotropic liquid the nearest neighbouring of a given molecule (the first 
coordination sphere) has a definite local symmetry, which can be characterized by 
the set of bond directions. The local structure of the liquid in the neighbourhood 
of a molecule at the point r' Q is characterized by the bond directions r' = r 2 — r' . It 
occurs that if the point r' Q is at sufficiently large distance from r then there is no 
correlation between the directions r = ri — ro and r' = r 2 — r' . In this case after the 
averaging over the system as a whole the pair distribution function transforms into 
the RDF and the equation ([[]) for s = 1 has the solution F 2 (ri — r ) = g(|ri — r |), 
which corresponds to ordinary isotropic liquid. When we approach the anisotropic 
liquid phase the long-ranged correlations between the bond directions r and r' do 
appear and the averaged two-particle distribution function depends on the bond 
direction now. 
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In the vicinity of the transition one can write 



F 2 (r 1; r ) = <?(|n - r |)(l + /(n - r )) (7) 

where /(n — tq) has the symmetry of the local neighbourhood of the particle at ro. 
The bifurcation point is given by the linearized equation ([[]) for s — 1, namely, 

/(r a - r ) = J r(r l5 r , r 2 )/(r 2 - r ) g{\r 2 - r \)dr 2 , (8) 



where 
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At the same time, when one approaches the line defined by the bifurcation con- 
dition, the correlation radius for the orientation fluctuations of the pair distribution 
function diverges. This fact can be shown with the use of the gradient expansion 
technique in the case of the equation ([]]) for s — 3, if we write the long range part 
of the correlator using the principle of vanishing correlations (Q) as: 

F 4 (r 1; r 4 ) = - r 2 |M|r 3 - r 4 |)(l + / 4 (n, .., r 4 )) (10) 

/ 4 (n, r 4 ) = / 4 (r, R, p, <p x , <p 2 ). 

Here (p± is the angle between the vector r = n ~~ r 2 an d the axis R = r 2 — r 3 , <p 2 is the 
angle between the vector p = r 3 — r 4 and the same axis. We have / 4 (r, R, p, (fi, <p 2 ) — > 
when R — > oo . 



The microscopic expressions for the elastic moduli and Frank constant [TI] 
enable us to understand on the microscopic level whether the 2D melting for any 
given potential is 3D like or whether it follows the KTHNY scenario. 

Let us consider now a possible description of the liquid-glass transition in terms 
of space symmetry breaking for three (four) particle distribution function in 2D (3D) 
systems. At high temperature the nearest neighbours of a molecule can take differ- 
ent relative positions and there is no short-range order (SRO). At lower temperature 
a SRO appears which can be of different kind at different densities (for phase transi- 



tions in liquids see ||11|| ). The rotation and the translation of the clusters of prefered 
symmetry give rise to the fact that one-particle and two-particle distribution func- 
tions remain isotropic. If a kind of bond orientational order (BOO) appears the 
clusters are oriented in similar way and the two-particle distribution function be- 
comes to be anisotropic (as in 2D hexatic phase). However, we can imagine another 
situation - freezing of the symmetry axes of the clusters in different position. The 
isotropic phase can be considered as analogous to the paramagnetic phase (of cluster 
symmetry axes), the BOO phase - to the ferromagnetic phase, and the mentioned 
freezed phase - to a spin glass phase. 
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Let us consider for simlicity a 2D system. In the vicinity of the transition one 
can write (in the superposition approximation for the liquid) 

F 8 (n|rS, rjj) = 0(|n - rS|)p(|ri - rg|)(l + / 8 (n|rS, rjj) (11) 

In 2D case /3(i*x|r°, r 2 ) depends in fact on two distances and two angles 

/ 3 (r 1 |r;,rO)=/ 3 (i2o,0o;i2i,e 1 ), (12) 

where Ro = r 2 — r°, Ri = T\ — r°, R2 = 1*2 — r? and 0o is the angle of the vector 
R with the z axis, 61 - the angle between Ri and Ro and O2 - the angle between 
R 2 and R . 

The linearization of (|I|) for s = 2 gives: 

/sORo, O ; 80 = | r'(i? , O ; r 2 ; e^/s^o, 0o; #2, e 2 )<?(|R 2 - K \)g(R 2 )dr 2 , 

(13) 

where 

r / (i2 O) 0o;r 2 ;-Ri,@i) = 51 ?! 7T7 / ^+1(1-1, r fc -4-i)5 , (|r 3 — r?|) 
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x ^(|r 3 - r°|)...s(|r fc+1 - r°|)<?(|r fc+1 - r°|) dr 3 ...dr k+1 (U) 

There are two kinds of angles entering the equations and two kinds of order pa- 
rameters, cosequently. One angle (0o) fixes the position of one pair of particles of 
the cluster, and the other (Oj) - the position of the third particle in the coordinate 
frame defined by 0o- The order parameter connected with 0j is the generalization of 
intracluster hexatic parameter for the case of different coordinate frames. The order 
parameter connected with O is an analogue of magnetic moment and in glass-like 
phase one can consider an Edwards- Anderson parameter < cos (f> (t) cos0 o (O) >• 
In such a way we come to the concept of a "conditional" long range order: if we 
consider two pairs of particles at infinite distance from one another then there exists 
a preferable possibility for the relative position of the third particle near each pair. 
The directions of the bonds in the pairs of particles themselves are subjects to spin- 
glass-like order. In 3D case the rotation of clusters is given by matrices D\™ (aJ i) 
so that we obtain a kind of orientational multipole glass for the clusters. 

This work was partially supported by Russian Foundation for Basic Researches 
(Grant No. 98-02-16805). 
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